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Abstract 

The Chazy equation y'" = 2yy" — 3y'^ is derived from the automorphic properties 
of Schwarz triangle functions S{a, (3,^; z). It is shown that solutions y which are 
analytic in the fundamental domain of these triangle functions, only correspond to 
certain values of a, (5, 7. The solutions are then systematically constructed. These 
analytic solutions provide all known and one new parametrization of the Eisenstein 
series P,Q,R introduced by Ramanujan in his modular theories of signature 2,3,4 
and 6. 



1. Introduction 



In a series of papers between 1909 - 1911, J. Cliazy considered a class of nonlinear 

differential equations of the form 

y'" - 2yy" + Sy'^ = ^^^(6y' - y')\ < M 6 , (1) 

during the course of his work on the extension of Painleve's program to equations of 
third order. The general solutions of these equations, which are now referred to as Chazy 
class XII, can be parametrized by solutions of appropriate hypergeometric equations for 
/c > 0, and the Airy equation when = 0. Furthermore, for 6 < G Z the solutions 
evolve from a generic initial data to form a natural boundary which is closed curve in 
the complex plane beyond which the solutions can not be analytically continued. For the 
special cases of = 2, 3, 4 and 5, the associated hypergeometric solutions are algebraic 
functions classified by Schwarz |4j, and leads a 3-parameter family of rational solutions 
for (pP). The focus of this note is on the equation 

y"'-2yy" + y = Q, (2) 

which corresponds to the limiting case ^ oo of ([1]), and will be referred to as the 
Chazy equation throughout this article. Chazy [2j showed that the general solution of 
([2]) also possesses a movable natural boundary by relating the solution y{z) to that of the 
hypergeometric equation s(s — l)x" + ~ ~ liiX = 0. 

In recent years, it has been shown that the Chazy equation arises in several areas of 
mathematical physics including magnetic monopoles [5], self-dual Yang- Mills and Einstein 
equations [HI El, and topological field theory [S]. In addition, ([2]) has been derived as spe- 
cial reductions of hydrodynamic type equations [H] as well as stationary, incompressible 
Prandtl boundary layer equations [TU] . These results have renewed interest in the study of 
the Chazy equation. For example, the SL2(C) symmetry of ([2]) was exploited to system- 
atically derive its general solution in [11]; moreover, in Ref. [12] the group invariance was 
applied to elucidate the role of its pole singularities, which under suitable perturbations 
coalesce to the natural boundary. Yet another interest in Chazy equation stems from its 
well-known connection with the automorphic forms associated with the modular group 
SL2(Z) and its subgroups. The notation r(l) = SL2(Z) for the full modular group (cf. 
[36] ) will be used throughout this article. It is quite significant that the emergence of the 
Chazy equation ([2]) in the theory of modular forms and elliptic functions can be traced 
back to the work of Ramanujan and subsequently in the work of Rankin and others. Here 
we briefly recall some of these interesting results which were anteceded by Chazy's work, 
but were apparently not noticed by those researchers. 

In 1916, Ramanujan [13], [HI pp 136-162], introduced the functions P{q), Q{q) and 
R{q) defined for |g| < 1 by 

oo„ 003^ 005^ 

n=l ^ n=l ^ n=l ^ 

(3) 

and proved by using trigonometric series identities that these functions in (|3]) satisfy the 
differential relations 

5P=^^, SQ=^^, SR=^^, S:=q-. (4) 
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By expressing the system (jl]) as a single differential equation for P{q), then defining 



y{z) := mP{q) , q ■= e^^" , Im(z) > , (5) 

one easily recovers the Chazy equation ([2]) [15j. (Note that in this case the natural 
boundary corresponds to the unit circle |g| = 1 or the real axis lm{z) = 0). Ramanu- 
jan's P{q), Q{q), R{q) correspond to the (first three) Eisenstein series associated with the 
modular group r(l), but this modern terminology is not needed in the present context. 
Ramanujan also considered the modular discriminant function 



n=l n=l 



he proved as well as conjectured many properties associated with the integer coefficients 
T{n) above, which are referred to as Ramanujan's tau-functions. It follows from the last 
two equations in (jlj) that P{q) is simply the logarithmic derivative of A(g), that is 

P(,) = ^. or yW = ^^. (6) 

Equation also follows by taking the logarithmic derivative of the infinite product 
formula for A(g) above, and comparing it with the g-expansion for P{q) in ([3]). Rankin 
in [Hj, showed using properties of modular forms that A{z) satisfies the equation 

2A""A^ - 10A"'A'A2 - 3A" ^ A^ + 24A"A' ^A - 13A' ^ = , 

which is homogeneous of degree 4 (in both A and its derivatives). Rankin's A-equation 
follows from the Chazy equation ([2]) with y{z) as in ([6]). In Ref. [16], Rankin also derives 
(among others) two more equations which are equivalent to 

AQ6^Q - b{5Qf = 960A , 6R6^R - 7{5Rf = -3024gA , 

which can be also deduced from the Ramanujan equations (jlj). Moreover, re-expressing 
the first equation above in terms of P{q) and its derivatives by employing (jlj) and the def- 
inition of A(g), yields (j2j) once more; while the second equation turns into the differential 
consequence of the Chazy equation. The first equation for Q above, was also obtained by 
B. van der Pol, who used it to derive the arithmetical identity [T7| 

n-1 

r(n) = n^cTa (n) + 60 ^^(2r2 — 3m) (n — 3m)o"3(m)o"3(n — m) , (1^(72) := c?'^ , riEN, 

m=l d\n 

relating Ramanujan's tau-functions and the sum-of-divisor function cr3(n). This and such 
other arithmetical identities follow from equating the g-expansions of both sides of a 
differential (or polynomial) relation involving the modular functions. We note here that 
using (jSj) in the Chazy equation leads to the identity 

n-1 

n'^{n — l)(Ji{n) + 12m(5m — 2n)ai{m)ai{n — m) = , 

m=l 
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that follows from the g-expansion of P{q) in ([3]). 

Ramanujan extensively studied the properties of his modular functions P,Q,R and 
established numerous identities involving them ITS] . Quite remarkably, he like Chazy, 
also employed the theory of hypergeometric functions to establish an implicit parametriza- 
tion of the functions P,Q,R, which plays a crucial role in the proofs of Ramanujan's mod- 
ular identities. In his second notebook [18], Ramanujan considered the hypergeometric 
function 2-^1(^5 1; 1; a^) related to the complete elliptic integral of the first kind by (see 

e.g., mm) 

K{x) := r ' , , = J 2F,a, i; 1; x) , 



JO — X sm t ^ 

and gave the following implicit parametrization of the functions P,Q,R 

P(q) = (l-5x)x^+12x(l-x)xx' , Q{q) = {l+Ux+x'^)x^ , R{q) = {l+x){l-34:X+x^)x^ . 

(7) 

where x{^) '■= 2-^1(5, h 1; and the nome q is defined by 



^2' 2 

2-^ 1 V 2 ' 2 



q = e , u := Tc 



2Fi{l,l;l;x) 



Equation ([7]) originates from Jacobi's work on elliptic functions ^19j. The derivation of ([7]) 
involves expressing P,Q,R as logarithmic derivatives of the quotients of theta functions 
and application of the remarkable Jacobi- Ramanujan inversion formula 



i?3(0|g):=l + 25^g"' = ^ 



_^x) 

n=l 



It is worth pointing out here that in Refs. [H [3], Chazy had also noted the equivalence 
between ([2]) and a system of three first order equations introduced by G. Darboux [20] in 
1878 (see Section 2). This first order system is different from the Ramanujan system (jlj) 
even though its solutions were given by G. Halphen [21] in terms of logarithmic derivatives 
of the null theta functions '^2,^3,^4. Each of the null theta functions satisfies Jacobi's 
nonlinear equation 



32^'" - i5ee'e" + 3oe'^y = 32(69" - 3^'')=^ + Tr^e'^^iee" - se'y = o , 



for 6{z) := '(9j(0|g) , 2 = 2... 4 and where q = e'^*^ p2]. Ramanujan relied heavily on the 
parametrization ([7]) and the inversion formula to develop his theory of modular equations 
involving functional relations among complete elliptic integrals at different arguments. 
Interestingly, he also proposed alternative parametrizations for P{q),Q{q), R{q) in terms 
of other hypergeometric functions (see Section 4), and where the appropriate nomes q can 
be expressed by 

a = e-"^ u ■= 2Fi(^,^;l;l-x) ^^ 2 3 4 6 (8) 

The index r is referred to as the signature of Ramanujan's theories. The case r = 2 
corresponds to Ramanujan's original theory of modular equations while signatures 3,4 and 
6 correspond to Ramanujan's alternative theories. He stated the results in his alternative 
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theories without proof [T8]. Proofs were constructed much later [23l [23] and is now 
an important topic of research [23]. Indeed a unified framework of Ramanujan's modular 
equations in different signatures based on the more contemporary theory of modular forms, 
elliptic surfaces and Gauss-Manin connections was proffered only recently in Ref |26] . 

In this paper, we derive the parametrization for the functions P,Q,R in Ramanujan's 
original and alternative theories, as well as some new ones (see cases (1), (2) in Table 
1) via the parametrization of the Chazy solution y{z) in terms of Schwarz triangle and 
hypergeometric functions. Our approach is not number theoretic, but rather based on 
the theory of Fuchsian differential equations and the action of its projective monodromy 
group on certain differential polynomials. Section 2 provides some necessary background 
on Fuchsian equations and their role in the conformal mappings of the complex plane 
to triangular domains T bounded by circular arcs. The triangle functions introduced 
by Schwarz form the natural coordinates on T, and which remain invariant under the 
group of automorphisms induced by the projective monodromy groups of the Fuchsian 
equations. The relationship between the triangle functions and the Chazy equation is 
established in Section 3. In section 4, we give the explicit parametrization of the Chazy 
solution y{z) in terms of the appropriate hypergeometric functions, thereby connecting 
our results to those of Ramanujan's parametrization in his theories of modular equations. 
Along these lines, we also note the work in Ref. ^27j which utilizes arguments motivated 
by the Lie symmetry analysis of (jlj) and related differential equations to obtain various 
parametrizations for the functions P,Q,R. The approach in this note is different from 
that in Ref. [27]. 

2. Conformal mapping and triangle functions 

The mapping properties of Fuchsian differential equations play a significant role in the 
theory of automorphic functions. In particular, H. A. Schwarz in 1873 carried out an 
exhaustive study of the conformal maps induced by ratio of solutions of hypergeomet- 
ric equations and the so called triangle functions [1]. This beautiful theory has since 
been developed significantly and is treated in numerous monographs. The brief overview 
presented in this section closely follow the texts [28] and [29] . 

A second order, Fuchsian differential equation with three regular singular points in 
the complex plane can be cast into the form 



where a, /3, 7 are the exponent differences (for any pair of linearly independent solutions) 
prescribed at the singular points 0, 1 and 00, respectively. (Note that /' indicates deriva- 
tion with respect to the argument of the function / throughout this article). The ratio 
z{s) of any two linearly independent solutions mi,M2 of ([9]) is a multi- valued function 
branched at the regular singular points, and satisfies the Schwarzian differential equation 

{z,s} = ^, {^^s}:=---^-j . (10) 

The monodromy group G C GL2(C) resulting from the analytic extensions of the pair 
(mi, U2) along all possible closed loops through an ordinary point Sq, is determined (mod- 
ulo conjugation) by the exponent differences. G acts projectively on the ratio z{s) via 
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fractional linear transformations 



z 7(z) : 



az 



cz + d 



7 = 



a b 
c d 



Since both 7 and A7 yield the same fractional linear transformation for any complex 
A 7^ 0, the projectivized monodromy group is the quotient group F = G/AJ C PSL2(C) 
where / is the 2x2 identity matrix. Note that both z and 7(2;) satisfy the same equation 
ffTU]) due to the invariance property of the Schwarzian derivatives: {z, s} = {7(2;), s}. A 
special class of solutions of flTU]) was extensively investigated by Schwarz who considered 
the parameters in V{s) to be real and < a,/?, 7 < 1. If V{s) is to be further restricted 
such that a + /3+7 < 1, then a branch of z{s) maps the upper-half s-plane (Im s > 0) onto 
a hyperbolic triangle T in the extended complex plane, bounded by three circular arcs 
which enclose interior angles avr, (3it, and ■jn at the vertices ^(0), ^(1), and z{oo). By the 
Schwarz reflection principle, the analytic extension of this branch to the lower-half plane 
across a line segment between any two branch points, maps the lower-half plane to an 
adjacent triangle T' that is the image of T under reflection across the circular arc which 
forms their common boundary. Continuing this process, the complete set of branches of 
z{s) maps the s-plane onto a Riemann surface spread over the z-plane consisting of an 
infinite number of circular triangles obtained by inversions across the boundaries of T and 
its images. The necessary and sufficient condition that this Riemann surface is a plane 
region D (being the uniform covering of non-overlapping triangles) is that the exponent 
differences a, j3, 7 be either zero or reciprocals of positive integers. In this case, the inverse 
s{z) is a single-valued, meromorphic, automorphic function whose automorphism group 
is the projective monodromy group V defined above. That is, 5(7(2;)) = s(z),7 G F. 
In this setting of conformal mapping, F is a discrete subgroup of PS'L(2,M), and turns 
out to be the group of fractional linear transformations generated by an even number of 
refiections across the boundaries of the circular triangles. More precisely, let ra,rf3,r^ be 
the refiections across the sides opposite to the vertices z{0), z{l), and z{oo) of T. Then the 
automorphism group is generated by the elements Ra = r^r^, Rp = Var^ and R^ = raVp 
which are rotations about each vertex by 27ra, 211(3 and 2717, and satisfy 



A vertex with a nonzero (interior) angle vr/m, m G Z"*" is called an elliptic fixed point of 
order m, whereas a vertex with zero angle is called a parabolic fixed point of the group. 

In its domain of existence D, the only possible singularities of s{z) and its derivatives 
are located only at the vertices where s{z) takes the value of 0, 1, or 00. Due to the 
automorphic property, it is sufficient to consider the function s{z) on the triangle T. Let 
2 = 2:0 be a vertex of the triangle T such that s{zq) = Sq G {0, 1, 00} is a regular singular 
point of the Fuchsian equation ^ with exponent difference G {«,/?, 7}. The behavior 
of the s{z) near a vertex zq depends on whether zq is an elliptic or a parabolic fixed point 



(i) If /i = l/m,m G then zq is an elliptic fixed point of order m. In this case, a 
pair of fundamental solutions of ([9]) in the neighborhood of s = sq 7^ 00 are of the form 



where ipi{s), i = 1,2 admit convergent power series in the neighborhood of s = sq, and 
'4'i{^o) 7^ 0. By taking appropriate linear combinations of the solutions ui and U2, z{s) is 



— R]/^ — RI^'^ — RaRfsRf — 1 • 



of F. 



u,{s) = {s - soY'+^'>/'Ms) , 



U2{S) = {S - SoY'-^'>^'Ms) , 
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defined via ^ 

where ipi^s) is analytic near s = sq and ip{s) ^ 0. The inverse function is single-valued, 
and given by 

s(^) = So + (^-^or0i(^), 01(^0)7^0, (11a) 

where (j){z) is analytic near z = zq. Thus, s — Sq has a zero of order m at z = Zq. If 
Sq = 00, then by making the transformation s' = 1 /s in ([9]) one finds in a similar manner 
as above that s{z) has a pole of order m at 2; = zq. That is, 

S{Z) = {Z- Zo)~™02(^) , M^O) ^ , (lib) 

and 02 (-2) is analytic in the neighborhood oi z = zq. (ii) When zq is a parabolic vertex, 
fi = 0. In this case, a pair of linearly independent solutions of ([2]) in a neighborhood of 
s = So 7^ 00 is given by 

Mi(s) = (s - soY^^tps^s) , U2{s) = [k\og{s - So) + ^4(s)]mi(s) , 

where '03(so) 7^ 0, '04(so) = 0, A; is a constant, and both ■03(s), '04(s) admit convergent 
power series in the neighborhood of s = So- In terms of Ui,U2, the function z{s) can be 
defined as 

^ . U2 2m U2 
Zmz = — or = — , 

Ml Z — Zq Ui 

depending on whether zq = 00, or a finite vertex in the extended ^-plane. From above, 
the inverse function s{z) can be expressed as a power series 

00 

siz) =so + J2 . ^ ■■= e'"*'^' (or <1 ■■= , c„ G C , (12a) 

n=l 

in the uniformizing variable q. Thus, s{z) is a single-valued function of z, holomorphic 
at g = 0. If So = 00, then a similar analysis as above can be carried out by making the 
transformation s' = 1/s, and by introducing the same local uniformizer q as in (I12al) . In 
this case, s{z) has a pole at g = 0, and the g-expansion 

, 00 

s{z) = - + J2dnq" , d,dneC. (12b) 

^ n=0 

A pair of adjacent triangles T and T' form the fundamental region X of the auto- 
morphism group r whose action on X tessellates the region D. The inverse function s{z), 
which maps X to the entire extended s-plane, generates the function field (over C) of 
X. The boundary of D in the 2;-plane, is a F-invariant circle which is orthogonal to (all 
three sides of) the triangle T and all its reflected images. This orthogonal circle is the 
set of limit points for the automorphic group F, it is a dense set of essential singularities 
forming a natural boundary for the function s{z). In its domain of existence D, the only 
possible singularities of s{z) are poles which correspond to the vertices where s{z) = 00. 
Thus, F is a Fuchsian group of the first kind (see e.g., [28], Sec. 30 for a definition) and 
s(z) is a simple automorphic function of F. Fuchsian groups associated with differential 
equations ([9]) with three regular singular points are referred to as triangle groups and 
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s{z) := S{a, P,'y; z) are called Scliwarz triangle functions. It follows from (fTOj) that the 
inverse function s{z) satisfies the following third order nonlinear equation 

{s,z}+'-^Vis)=0. (13) 

Conversely, when the parameters 7 in V{s) are either zero or reciprocals of positive 
integers, a three-parameter family of solution of f|T3l) is obtained as the inverse of the 

pK±Bu^ AB.C.DeC. AD-BC=1, (14) 

Cmi(s) + Du2[s) 

where Ui and U2 are linearly independent solutions of ([9]). The solution is single- valued and 
meromorphic inside a disk in the extended 2;-plane, and can not be continued analytically 
across the boundary of the disk. This boundary is movable as its center and radius are 
completely determined by the initial conditions which depend on the complex parameters 
A,B,C,D. 

A number of nonlinear differential equations whose solutions possess movable natural 
boundaries, can be solved by first transforming them into a Schwarzian equation ( |T3i) 
and then following the linearization scheme described above. We briefly recount some 
examples of such nonlinear differential equations and their relations to Schwarz triangle 
functions. In 1881, Halphen considered a slightly different version [30l pp 1405, Eq (5)] 
of the following nonlinear differential system 

w'l = -W2WS + Wi{w2 + ws) + , 

W2 = -W3W1 + W2{W3 + Wi) + t'^ , (15) 
W3 = -W1W2 + ws{wi + W2) + , 

= a^{wi - W2){W2 - W3) + 0^{W2 - Wi){wi - W3) + 7^(W3 - Wi){w2 - W3) , 

for functions Wi{z) 7^ Wj{z), i 7^ j, i,j G {1,2,3}, and constants a, P, 7. Halphen pre- 
sented the solutions of (1151) in terms of hypergeometric functions. More recently, the 
authors found that (|T5|) arises as a symmetry reduction of self-dual Yang-Mills equa- 
tions, and called it the generalized Darboux-Halphen (gDH) system [311 [32]. If in fact, 
Wi{z), W2{z), W3{z) are parametrized in terms of a single function s{z) (and its deriva- 
tives) as 



ir, /s 



log 

s 



1 

W2 = - 



log ' 



s-1 



1 

W3=- 



log 



s(s-l] 



(16) 



then s{z) is a solution of ffT^ . where the constants a, P, 7 in V{s) of ffTSl) are the same 
as those appearing in of ffT^ . The special case a = /5 = 7 = 0of the gDH system 
corresponds to the "classical" Darboux-Halphen (DH) system which is equation (fT5|) with 
= 0. This equation originally appeared in Darboux's work of triply orthogonal surfaces 
on in 1878 [20], and its solution was subsequently given by Halphen [21] in 1881. 
The variables Wi associated with the DH system can be parametrized as in (1161) by the 
triangle function s{z) = S{0, 0, 0; z); the latter is related to the elliptic modular function 
X{z) := -(92 (Ok) 7^3(0^)5 expressed in terms of null theta functions. Chazy showed that 
the function y{z) := 2{wi + W2 + W3) satisfies ([2]) introduced in Section 1. Chazy [3] 
also noted that besides the elhptic modular function X{z), the solution to ([2]) can also 
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be given in terms of the triangle function S'(|, |, 0; 2;) which is the same as the modular 
J-function for the group SL2(Z). In more recent work Bureau [33] in 1987, investigated 
a class of third order nonlinear equations, and expressed their general solutions in terms 
of the Schwarz triangle functions s := ^(a,/?, 7; z). In particular. Bureau's class includes 
the Chazy equation ([2]). 

It is natural to inquire whether the solution of the Chazy equation ([2]) admits parametriza- 
tion in terms of ot/ier triangle functions besides S{0,0,0]z) and S{^, ^,0; z). One mo- 
tivation of the present work is to address this question and to investigate the possible 
linearizations of the Chazy equation ([2]) via solutions of the Fuchsian equation (Q with 
parameters {a, P, 7} other than {i, |, 0} or {0, 0, 0}. In the following we outline a method 
to systematically derive the Chazy equation from the solution of the Schwarz equation 
( IT3l) for appropriate values of the parameters a,/3,7. Our construction utilizes the trans- 
formation properties of s{z) and its derivatives under the automorphism group F and 
selects those triangle functions which provide a natural parametrization of the Chazy 
solution y{z) that is holomorphic in their domain of existence D. 

Triangle functions and the Chazy equation 

Let F C PSL2(M) be a Fuchsian triangle group with a,/3, 7 either zero or a reciprocal of 
positive integers, as in the previous section, and let s{z) be a simple automorphic function 
of F defined on a domain D of the complex plane. A meromorphic function / on D is 
called a automorphic form of weight k for F if 

/(7(^)) = {cz + dff{z) , 7 = ^) e F , 7(^) = ^ 

for all z G D. If A; = 0, then / is called a automorphic function on F and is a rational 
function of s{z). Consider a gDH system (fT5l) for F where the gDH variables Wj, i = 1, 2, 3 
are parametrized as in f|T6l) . From the automorphic property: 5(7(2;)) = s{z), it follows 
that ^'(7(2)) = {cz + (i)^s'(z), and that 

Wi{'y{z)) = {cz + dfwi{z) + c{cz + d) , 7 G F . 

That is, s'{z) is a weight 2 automorphic form, whereas the Wi are called gwasz-automorphic 
forms of weight 2. Define in terms of the gDH variables, the following function 

y{z) = aiwi + a2W2 + a^Ws , (17) 

on D, where the coefficients are constants. The objective of this section is to determine 
an autonomous differential equation which is a polynomial in y and its derivatives. 

It follows from the transformation property of the Wi above, that y{z) transforms 
under the action of F as 

y{-f{z)) = {cz + dYy{z) + pc{cz + d) , 7GF, (18) 

where p = ai -|- 02 + 03 is called the coefficient of affinity of the quasi-automorphic form 
y{z). Furthermore, a sequence of automorphic forms can be constructed on the ring of 
differential polynomials of y{z) as follows: 
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Lemma 1. Let y{z) be a quasi- automorphic form of T with affinity coefficient p, then 
f2 = y' — y'^/p is an automorphic form of weight 4, and fn+i = f'n~ {'^''^/p)yfn, n>2 are 
automorphic forms of weight 2n + 2 ofT. 

Proof. Differentiating f|T8|) . one finds that 

y'i^iz)) = {cz+df {{cz + dfy{z) +pc{cz + d))' = {cz+df [y'{z) + |^ + 



and also from f|T8|) . 



cz ^ d [cz + d)" 

Tlien, by combining tlie two expressions above yields the desired transformation property 
for /2. The rest follows by induction and use of (JTSll . □ 



The parametrization f[T6|) of the gDH variables and (jTTjl yield the following expression 
for yi^z) in terms of and its derivatives 

yiz) = ^-Z^ Mz) = - 19 

' 2 0(2) ' ' {s- 1)^1 s''^ ' ^ ' 

with h\ = 1 — ai/p and 62 = 1 — 02/^. Then all the higher derivatives of y can be 

expressed also in terms of s{z),s'{z) and s"{z) by differentiating f|T9|) successively and 

using the Schwarz equation f|T3|) for s{z). In particular, the differential polynomials fk 
introduced in Lemma [1] are given by 

fn = ~Vn{s)s'{zr, Vn+i{s)=V;{s)+nq{s)Vn{s), n>2, (20) 
where Vn{s) are rational functions of s defined recursively from 

V^is) = '^V{s) + q'{s) + Iqisf , = + ^ , (21) 

^ ^ s — 1 s 

and where V{s) is given in (Q. Combining Lemma [1] with (!20!) for n = 2,3,4, and with 
/2 7^ 0, one can construct the following rational expressions in y{z),y'{z),y"{z) and y"'{z) 



fi 2Vs{sy /4 2^4(5) 



fi pV2{sY' fi pv^is) 



(22) 



which are rational in s{z). Eliminating s from the two equations in ( l22l) leads, in the 
general case, to a third order, nonlinear equation, rational in y{z) and its derivatives, and 
which depends explicitly on the parameters p, 6i,&2 (equivalently, 01,02,03), and a,P,'j. 
For suitable choices of these parameters, the third order equation constructed this way can 
be identified with a large set of nonlinear differential equations, and whose solutions y{z) 
are given via (1191) . in terms of the solutions of the Schwarzian equation (1131) . Furthermore, 
it follows from Section 2 that all such solutions will admit a natural boundary, and will 
be meromorphic in the domain of existence D. Of these third order equations, only the 
special case of the Chazy equation will be considered in this note, leaving the general 
classification problem for a future work. 
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The Chazy equation ([2]) can be expressed simply as a polynomial equation in terms 
of the automorphic forms /2 and /4. Indeed, from Lemma [T] 




Then imposing the constraint 



p := ai + 02 + as = 6 , 



on the coefficients Oj in (ITTl) in the above expression for f^, yields the alternative expres- 
sion: /4 + 4/1 = for the Chazy equation ([2]). Note that this expression is equivalent 
to the vanishing of a certain automorphic form of weight 8 associated with the Fuchsian 
group r (see Section 4). Then the second equation in fl22|) implies that the Chazy equation 
must be equivalent to the condition 



which needs to hold for all s. Condition fl23|) imposes certain restrictions on the parame- 
ters 61,62, and a,/?, 7 appearing in V2(s),V4(s). However, we impose further constraints 
on the parameters 61,62 by demanding that the meromorphic function y{z) be in fact, 
holomorphic in its domain of existence D. The reason for this additional condition is mo- 
tivated from the known result (see e.g. [2]) that the general solution y{z) of the Chazy 
equation obtained via the triangle function S{^,^,0; z) is analytic on D although there 
exists particular (2-parameter family) solutions that are meromorphic but do not possess a 
natural boundary. Since one of the main objectives of this paper is to obtain parametriza- 
tions of the Chazy equation in terms of triangle functions, we impose the holomorphicity 
of y{z) a priori; then this leads to the specific choices for the parameters {a, (3, 7} for the 
triangle functions which parametrize the Chazy solution y{z). These algebraic conditions 
will be systematically investigated next. 

In order to determine whether y{z) defined in ( ITTl) is holomorphic on D, it is neces- 
sary to analyze the singularities of the gDH variables Wi given in terms of s{z) and its 
derivatives in (fT6l) . It follows from the conformal mapping theory discussed in Section 2 
that it is sufficient to examine the behavior of the function s{z) and its derivatives near 
the vertices ^(0),2;(1), and z{oo) of the fundamental triangle T. The Schwarz reflection 
principle and the automorphic property then ensure that s{z) will have the same behavior 
at the vertices of the reflected triangles in D. 

Lemma 2. Let s{z) = S{a,(3,'^\ z) he the Schwarz triangle function of a Fuchsian group 
r with fundamental triangle T whose interior angles at the vertices z{{]) ^ z{l) , z{oo) are 
respectively, an, [371 and 77?. Furthermore, let the gDH variables wi,W2,W3 be given by 
s{z) and its derivatives as in f|T6|) . //{a,/3,7} = {1/pi, l/p2, l/ps}, Pj G Z"*", j = 1,2,3, 
then wi,W2,W3 have first order poles at each of the vertices with the following residue 
scheme: 



ReS^=2(0){w^l,W2,W3} = {-|, |(Pl - 1), -|} , ReS^=^(i){wi,W2,W3} = {|(P2 - 1), -|, - 

^eS^=,(oo){wi, W2, W3} = {-i, -i, i(p3 - 1)}. 



V4 = 12Vi , 



(23) 



Proof. Direct computation of ^1,^2,^3 in (fT6|) using equations (lllal) and (lllbO . 



□ 
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It follows from Lemma [2] that only possible singularities of y{z) defined by (fT7|) in the 
domain D, are first order poles at the vertices corresponding to the elliptic fixed points of 
r. (Note however that y{z) has a dense set essential singularities at the boundary of D as 
is the case for the function s{z)). Therefore, y{z) will be holomorphic in the entire domain 
D if the coefficients a, in f|T7|) can be chosen such that the residue at each pole vanishes. 
Consider first the case {a,/3,7} = l/p2, l/Ps}, where pj G Z+. From part (i) of 

Lemma m the condition that the residue of y{z) vanishes at each pole z{0), z{l), z{oo), 
leads to a set of homogeneous, linear equations for a^, namely, 

ai + (1 - pi)a2 + as = , 

(1 -P2)ai + a2 + a3 = 0, (24) 
fli + ^2 + (1-^3)^3 = 0. 

For nontrivial solutions Oj, the vanishing of the determinant of the coefficient matrix in 
fl24|) implies that 1 /pi + l/p2 + l /ps = a + /5 + 7 = l, while the conformal mapping of the 
upper-half s-plane onto the hyperbolic triangle requires that a + /3 + 7 < 1. Therefore, 
only two out of three equations in ( 124|) can be used to set the residue equal to zero. Thus, 
for y{z) to be holomorphic on D, one of the three vertices must be a parabolic vertex at 
the boundary of D instead of an elliptic vertex with a pole singularity. Hence, we have 
the following result. 

Proposition 3. If a solution y{z) of the Chazy equation ([2]) is parametrized by the triangle 
function s{z) = S{a, z) defined on a domain D, then a necessary condition for y{z) 
to be holomorphic on D is that at least one of the parameters a, P, 7 be zero. 

Here we emphasize that it is indeed possible for the solutions y{z) of other nonlinear 
equations obtained from the general construction outlined below fl22|) . to be meromorphic 
in its domain of existence D. For example, the general solutions of ([1]) for 6 < G Z 
are given in terms of the triangle functions 5(1, |, 2;), and are meromorphic in D with 
poles at the elliptic vertices which have interior angles 7r/k [2]. We will consider the 
parametrization of ([1]) and other equations in a future study. 

According to Proposition [3l at least one of the vertices of the fundamental triangle 
T must be a parabolic fixed point of the automorphism group F. Thus there are three 
distinct cases which are considered below. For reasons that will be clear in Section 4, z{0) 
is always chosen to be a parabolic vertex in each of these cases so that the corresponding 
triangle functions are of the form S{0, f3, 7; z). 

2 elliptic and 1 parabolic vertices: As indicated above, it suffices to choose z{fS) 
to be the parabolic vertex, and z{l),z{oo) as the elliptic vertices so that {a,/3, 7} = 
{0, l/p2, l/Ps}- All other sub cases can be generated from permutations of {-2(0), -2(1), z{oo)} 
Requiring the residue of y{z) to be zero at each of the poles -2(1), ,2(00) yields the 
last two equations of From these, one easily deduces that Oi = p/p2 = Pp and 

% = p/P3 = IP so that a2/p = 1 — — 7. Moreover, the parameters 61, &2 in f|T9|) are 
given by 61 = 1 — ai/p = 1 — /5 and 62 = 1 — ^2/^ = P + 1- Equations (|2T|) and fl20l) then 
deliver the explicit forms for V2(s) and V4(s), which depend on the remaining parameters 
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/?, 7- They are given by 

V2 = -^^y With A = and, 

_ (1 - 27)(1 - 37)g^ + [{1-13- 7)(127 - 7) + (1 - 27)]^ + 6(1 - - 7)^ 
^ s\s-lf 

Finally, inserting the above expressions in fl23|) and requiring that the resulting expression 

s[(l - 27)(1 - 87)5 + (1 - /5 - 7)(67 - 6/3 - 1) + (1 - 27)] = 

holds for all s ^ {0, 1, 00}, give rise to two distinct solutions: 

{a, /?, 7} = {0, i, |} , {a, A 7} = {0, |, |} . 

1 elliptic and 2 parabolic vertices: Again, without any loss of generality, the elliptic 
vertex can be chosen as z{l) such that {a,/?, 7} = {0,l/p25 0}. From Lemma [H y{z) 
is holomorphic in T except at the vertex z{l) where it has a simple pole. Vanishing 
of the residue at the pole z{l) gives the second equation in (12^ . which implies that 
Oi = p/p2 = Pp and bi = 1 — ai/p = 1 — p. In this case, 62 and (3 are the two remaining 
parameters in V2(s) and V4(s), which take the forms 

y _ {l-b2? B _ N{s) 

' 2s2 ^s(s-l)' ^ s4(s-l)3' 

where B := {1 — l3){2b2 — (3 — l)/2 and N{s) is a cubic polynomial whose coefficients 
depend on /5, 62- Then fl23l) implies that Cq — Cis{s — 1) = for all s ^ {0, 1, 00}, where 

Co = (i^MizMii^, c, = (l-/3)(2-3/?)(l-62)^ + S[6(l + /32)-5(62 + /5)] , 
and B is defined above. Setting Cq = = Ci, three distinct solutions are found: 

(0{«,A7} = {0,i,0}, 62G{i,|}; (M){a,A7} = {0,|,0}, 62 e {i, §}; 

(m){a,/?,7} = {0,|,0}, 62 = !. 

Notice that in case (in), /3 is not reciprocal of a positive integer, hence 5(0, |, 0; z) is not 
a single valued function of z. 

3 parabolic vertices: Here, 7} = {0,0,0} so that y{z) is holomorphic on D for 

any choice of the parameters a,. The functions V2{s) and V4(s) depend on the parameters 
61, 62 which are determined from fl2^ by requiring that 7/(2;) satisfies the Chazy equation 
(12]). The calculations are tedious but similar to the previous two cases, and there are two 
distinct solutions: &i = &2 = |, and 61 = 62 = §• 

The triangle functions S{a, (3, 7; z) corresponding to the three cases considered above 
are the ones which give holomorphic solutions of the Chazy equation ([2]) via (fT6|) and (ITTj) . 
These results are summarized in Table 1 below. Recall that 01 + 02 + 03 := p = 6, and 
that y = 3(j)' (z) / (j){z) as given in (1191) . The automorphic groups F corresponding to the 
triangle functions are listed in the second column of the table. With the exception of Case 
5, all others are subgroups of the modular group F(l) (see e.g., [36] for notation); these 
will be discussed in the following section. The triangle function in Case 5 is not a simple 
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automorphic function of a Fuchsian group of first kind since the exponent difference at 
the vertex z{0) given by /? = | is not reciprocal of a positive integer. Hence S{0, |, 0; 2;) 
is not a single valued function of z on the domain D. Nevertheless, one can show from 
the conformal mapping properties of the triangle functions that the single- valued function 
S{0, 5, z) can be expressed as a degree-2 rational function of 5(0, |, 0; z) [29], namely 



S{0,ll;ez) 



-45(0,1,0;;.) 
[5(0,|,0;^)-1]^ 



with e = ^ — |. Then (j){z) in case 5 is a constant multiple of the (f){ez) in case 1; hence, 
y(z) obtained from 5(0, |, 0; z) is a single-valued function. 



Case 


r 


s{z) = S{a,(3,-i] z) 


?/ = aiWi + a2U72 + ci^wz 


(piz) = s'{z)/{s-lf's^^ 


1 


r(i) 


'S'(0, |, -^'^ z) 


y = 3wi + W2 + 2w3 


= S'/{S - l)V2s5/6 


2 


r2 


*S'(0, 3, ^; z) 


y = 2wi + 2w2 + 


= S'/{S - 1)2/3^2/3 


3 


ro(2) 


5(0, 1,0;;.) 


(i) y = 3wi + 2w2 + ^i'3 

(ii) y = 3wi + W2 + 2w3 


(i) = s'/{s - l)V2s2/3 

(ii) = S'/{S - l)V2s5/6 


4 


ro(3) 


5(0, |,0;^) 


(i) y = 2wi + 3^2 + 

(ii) y = 2wi +W2 + 3ws 


(i) = s'/{s - ly/h^/'' 

(ii) = S'/{S - l)2/3s5/6 


5 




5(0,1,0;^) 


y = Awi + W2 + W3 


= s7(s-l)^/^s^/« 


6a 


ro(4) 


5(0,0,0;z) 


y = Wi+W2 + 4w3 


= S7(S - l)5/S5/6 


6b 


r(2) 


5(0,0,0;^/2) 


y = 2{wi + W2 + W3) 


= S7(S- 1)2/3^2/3 



Table 1: Triangle functions associated with the Chazy solution y{z) 

Note that in each of the cases 3 and 4, the two different forms of y{z) can be trans- 
formed to each other by interchanging W2 and W3. This transformation stems from the 
permutation of the two parabolic vertices z{0) and z{oo) of the fundamental triangle T, 
thereby inducing the inversion map s — » of the function field s{z). Under this involu- 
tion, 0(z) corresponding to the sub cases (i) and (ii) are constant multiple of each other, 
for both cases 3 and 4. In case 6, there are two automorphism groups which are conjugate 
to each other: ro(4) = gT{2)g~^ where g = (02)5 fundamental domain of r(2) is 

mapped under z ^ 2z to the fundamental domain of ro(4). The canonical automorphic 
function for r(2) is the elliptic modular function X{z), which is related to the triangle 
function of ro(4) as 5(0,0,0;z) = \{2z)[\{2z) — 1]^^ From the above relation, it turns 
out that the (f){z) in case 6b can be transformed to that in case 6a due to the well-known 
[see e.g., [29j) functional relation: X{z) = A^X{2z)[l + ^X{2z)]-^. Thus we see that 
there are many different representations of the solution of the Chazy equation in terms 
of Schwarzian triangle functions which in turn satisfy Schwarzian equations (1131) . As we 
have seen, each of these equations can be linearized via (Q. We turn to this topic next. 

4. Parametrization of the Chazy solutions 

In the previous section, explicit solutions of the Chazy equation were presented in terms of 
the triangle functions listed in Table 1. The standard expressions of the triangle functions 
are usually via the theta or Dedekind's eta functions which admit Fourier or g-series 
expansions (as in (I12al) or (I12bp ) developed in the neighborhood of the parabolic vertex 
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Zo = ioo (see e.g., [26] and references therein). It is however more convenient to express 
the Chazy solution y{s{z)) implicitly, that is, in terms of the variable s and a solution 
u{s) of the linear equation Q. It is then possible to treat the nonlinear Chazy equation 
purely on the basis of the classical theory of linear Fuchsian differential equations. This is 
the main purpose of the present section. Yet another motivation is to relate our results to 
Ramanujan's work on the various parametrization of his functions P, Q, R, as mentioned 
in the introduction. 

In this section, the domain of the automorphic functions s{z) will be taken as D = 7i, 
the upper-half complex plane, and the hypergeometric form of the Fuchsian differential 
equation ([9]) will be considered, in order to make contact with standard literature. If u{s) 
is a solution of ([9]), then the function 

;^(S) =s("-^)/2(s-1)(^-^)/2m(s) (25) 

satisfies the hypergeometric equation 

which can be written in more standard form as 

s{s - l)x" +[{a + b+ l)s - c]x' + abx = 0, (26b) 



where a = (1 — « — /5 — 7)/2, b = (1 — a — /5 + 7)/2, and c = 1 — a. The transformation ( l25i) 
sets one of the exponents to at each of the singular points s = 0, 1, but the exponent 
differences as well as the ratio z{s) of any two linearly independent solutions, remain 
the same as those in (Q. Consequently, the conformal mapping and the construction of 
the triangle function described in Section 2 can be carried out in an identical manner 
by employing the classical theory of the hypergeometric equation instead of ([9]). Let 
{Xi)X2} be a pair of linearly independent solutions of fl26ap or ( I26bl) . and set z{s) = 
X2{s)/xi{s), then s'{z) = l/z'{s) = xl/W{xi,X2) where W{xuX2) = Cs"-\s - 
is the Wronskian, and C 7^ is a constant depending on the chosen pair of solutions 
{Xi)X2}- Use of this expression for s'{z) in f|T9l) provides an implicit parametrization for 
y{z) in terms of Xi aiid its s-derivative, given by 

yi^s)) = |.-"(. - ir^{2s{s - l)xix[ - lih + h)s - h]xl) , (27) 

with bi = bi + (3 — 1 and 62 = ^2 + « — 1- In the following, we construct the triangle 
functions s{z) for all the cases listed in Table 1 as well as the corresponding Chazy solution 
y{z) using fl7r|) . In each case, the pair of hypergeometric solutions {xi, X2} are so chosen 
that the conformal map results in a fundamental region T which has the parabolic vertex 
z{0) = ioo G TC, around which suitable g-expansions for s{z) are developed. In addition, 
one imposes boundary conditions on y, y', y" in a consistent manner at this vertex in order 
to uniquely fix a special solution y{z), which is then shown to be related to Ramanujan's 
P{q) via The ensuing results then relate to Ramanujan's parametrization for his 
theories of signatures 2,3,4, and 6, as mentioned in Section 1. To the best of our knowledge 
the cases from Table 1 that were not recorded by Ramanujan, correspond to Case 1 which 
was known by Chazy himself [21 [3], and Case 2. The representation of y{z) via the 
triangle function ^(O, |, |; z) in case 2 can be found in Ref. 
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4.1. Case 1 



Here a = 0, /3 = | and 7 = |- The corresponding group of automorphisms, F, is generated 
by rotations about the vertices z{l) and z{oo) by tt and 27r/3, respectively, and a parabohc 
transformation stabihzing the vertex z{Q). It is known (see e.g. [28]) that the projective 
action of F is isomorphic to that of the modular group F(l) := SL2(Z) acting on the 
upper-half plane Ti via fractional linear transformations according to 

and which is generated by the fundamental transformations: 

z ^ z + 1 (translation) , z — (inversion) . 

z 

It is also customary to choose the fundamental triangle T in 7i such that it is a strip par- 
allel to the imaginary axis, bounded by the unit circle \z\ = 1, and the vertices are located 
at z(0) = ioo, z{l) = i, z{oo) = e^'^'^l'^ . To construct the triangle function S'(0, |, |; z) on 
T, one starts with the hypergeometric equation (cf. fl26ap . (126bp ) 

X"+ f- + ;T^^')x' + 4^^X = 0, (28) 
\s 2(s — 1)/ s\s — 1) 

with parameters a = ^, 6 = ^ and c = 1. Equation (!28|) admits a one- dimensional space 
of single- valued solutions spanned by the hypergeometric function 2-^1 (^, Ij -s) that is 
holomorphic in a neighborhood of s = and is normalized to unity there. The second 
independent solution of (125]) is chosen as 2-Fi(^, |i 1 ~ -s), which is holomorphic in a 
neighborhood of s = 1. Next, define the conformal map z(s) as follows 

z{s) = ^ , xi = 2Fi(^, h 1; s) , X2 = A 2Fi(^, A; 1; 1 _ s) + Bxx , (29) 
Xi 

where A, B are constants to be determined by fixing the vertices of T as specified above. 
The analytic continuation of 2-^1 (^, ^\\\^ ~ iiito the neighborhood of s = is given 
by (see e.g. [[35]), 



Ffi 

J-'^lv 19 ' 19 ' 9 ' 



12.12.2. . f(3^)f(A; 



n=0 



l0g(.)2F,(^,A;l;,) + 5^ii^/,„, 



where F(-) are Gamma functions, /i„ := 2-?/'(l + ra) — V"!]^ + Ti) — "^(^ + '^)) and ^/'(■) : = 
F'(-)/F(-) are the Digamma functions. Hence from (129|) . 



^(s) = A '^^^^V/'j'! + 5 = Ai(- log(s) + h{s)) + B, A, = A . ^^^^ 



^Fiihh^-.s) n v;; , F(^)F(A) ' 



near s = 0, and h{s) is holomorphic in that neighborhood. The constant A is determined 
by fixing Ai = i/2n above, so that z{s) — > ioo as s — > O"*", and z ^ z + 1 onto the next 
branch as s makes a circuit around s = 0. The constant B can then be determined by 
demanding that limz(s) = i. This sets the elliptic vertex z{l) = i, and yields B = —i 

s— >1 
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from the expression of z{s) above, after using 2Fi{j^, 1; 1) = r(|)[r(^)r(Y|)]~^ and 
the formula r(x)r(l — x) = 7rcsc('7rx). Consequently, the function q := e^'^" has the 
power series representation q = se-'i(«)+2T = Bis{l + ais + a2S^ + ...), which can be 
inverted to obtain a g-series representation of the triangle function s{z) = 5(0, |, |; z) in 
the neighborhood of the parabolic vertex q = 0{z = ioo). In particular, the constant Bi = 
g27r-/i.(o) evaluated from the the analytic continuation formula of 2Fi{j^, ^; |; 1 — s) 

given above. One easily checks that h{0) = Hq := 2?/'(l) — ip{j2) ~ ^(n) which evaluates 
to ho = 27T + 3 log 12 using Gauss's Digamma formula (see ^35]), and thus Bi = 1/1728. 
The first few terms of the g-expansion of the triangle function is then given by 

^(0, 1, 1; z) = 1728g(l - 7Uq + 356652g2 + ...). 

Using the analytic continuations of the pair {xi,X2} to a neighborhood of s = oo, and 
from the values of the constants A, B obtained above, it can be verified that z{oo) = e^'^*/^. 
We note here that S{0, z) is the reciprocal of the well known J invariant associated 
with the modular group r(l). 

In order to derive the implicit parametrization y{z{s)) for the Chazy solution, one 
needs to calculate the Wronskian W{xiy X2) of the two solutions specified in fl29|) . of the 
hypergeometric equation fl28l) . A short calculation employing the analytic continuation 
formula for 2-Fi(^, ^; |; 1 — s) near s = gives 

W{xi, X2) = ^[W{xu h{s)) - ^] = Cs-\s - \)-\ , 

where hi^s) is analytic near s = 0, and the last expression on the right follows from Abel's 
formula. Hence, letting s — > 0"^ in above, yields the constant C = l/27c. Then from fl271) 
with &i = I and &2 = | (cf. Table 1), the following parametrization is obtained 

vizis)) = mil - s)kxl + 12sxix'i) , Xii-s) = 2^1(3^, ^; 1; s) . (30) 

Recall from Section 1 that Ramanujan's modular function P(g) introduced in ([3]) is related 
to the Chazy solution via yiz) = niPiq). Therefore, from (l30l) it is possible to obtain 
an implicit parametrization for Ramanujan's P,Q,R in terms of the solutions of the 
hypergeometric equation (125]) . 

Proposition 4. Let zis) be the quotient of hypergeometric solutions xi o,nd xi defined in 
fl29l). andq = e'^^"^'\ then 

Piq) = (1 - s)kxl + i2sxix'i) , Qiq) = xt , Riq) = (i - s)\xl ■ 

Proof. First, from ([2]) note that the Ramanujan functions satisfy the conditions P 
1,(5— >1,-R^lasg— >0. Next, from ([2]) and using the forms /2, /s from Lemma [T] it is 
easy to verify that the triple {y/vr, ■^f2, j^^fz} satisfies Ramanujan's differential system 
dl]). On the other hand, it follows from (!30|) that y/vri^lass— >-0 (equivalently, q 
0). Moreover, differentiating yiz) in (!30!) successively and using s'iz) = Xi/^(Xi?X2), 
one obtains the expressions -^f2 = Xi and (-fya/a = \/l — s Xi, which satisfy the same 
conditions as Q and R, respectively when g — > 0. Therefore, uniqueness of solutions of 
the differential system (jlj) yields the desired result. □ 
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Substituting the g-expansion for S{0, ^, z) into the hypergeometric series for xi, it is 
possible to recover from the above parametrizations the g-series for P,Q,R in On the 
other hand, Proposition H] provides an elegant representation for the the triangle function 
S{0, 2, 3; -z) as well ClS db remarkable identity, namely 

q3 _ ^2 I 3 2 

5(0, \; z) = — — — , Q'i = 2Fi{jK, -j^; 1; ^753^) • 



The first expression can be used together with ([3]) to obtain the g-expansion for S{0, 2)3!^) 
derived above. 

The transformation property flTSl) for y{z) under the action of the modular group 
r = r(l) implies that P = y/iriis & quasi-modular form of weight 2 and affinity coefficient 
j9 = 6 on r(l). Moreover, by comparing (jlj) with the result of Lemma [1], it follows that 
Q = and R = j:^fz are modular ioims of weight 4 and 6 respectively, on r(l). In 
fact, P = E2, Q = E4 and R = Eq, where Ek{q) is the normalized Eisenstein series of 
weight k for the modular group r(l). E^ is a holomorphic modular form, and is defined 
(for even positive integer k) as 

2k ^ 

E,iq) = l- — J2^k-iin)q\ q = e'-^\ zeH, (31) 

where Bk is the kth Bernoulli number and (Tk{n) is the sum-of- divisor function introduced 
in Section 1. Detailed discussions of modular forms appear in several monographs (see 
e.g. [36l [37j. The vector space of holomorphic modular forms of weight k on r(l) is 
denoted by Mfc(r(l)). Proposition H] provides a parametrization for any / G Mfc(r(l)) 
since it is known from the theory of modular forms that / belongs to the polynomial ring 
C[E4, Eq\. For 4 < A; < 10, Mfc(r(l)) is one-dimensional, and is spanned by the Eisenstein 
series Ek- It turns out that the left hand side of ([2]) can be expressed as + 4/| (see 
Section 3) which by virtue of Lemma [H must be in M8(r(l)). Hence, /4 + 4/| = CEg. 
But with y = 7riE2, the expression /4 + 4/I is a differential polynomial in E2, and vanishes 
as g — *^ 0, whereas Eg 1. Therefore, the Chazy equation + 4/| = follows from this 
"modular" argument. 

4.2. Case 2 

The triangle group associated with a; = 0,/3=|,7 = |is denoted by F^. It is a normal 
subgroup of the modular group F(l) of index 2, generated by the period-3 transformations 

Ti : z^ —— , T2: , 

z + 1 z 

whose respective fixed points e^'^^^^ and e'"^^^ in Ti. form the vertices z(l) and z(oo) of 
the fundamental triangle T of F^ together with the parabolic fixed point z{0) = ioo. An 
important point to note here is that the stabilizer of the vertex z{0) = ioo in F^ is given 
by the transformation T2T1 : z ^ z + 2 instead of the translation z ^ z + 1, which is not 
an element of F^. 

One proceeds with the construction of the triangle function S{0, |, |; z) in a similar 
manner as in Case 1. Define the conformal mapping by 

zis) = ^ , xi= 2F,il i; 1; s) , X2 = A 2Fi(|, 1; |; 1 - s) + Bxi , (32) 
Xi 
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where 2Fi(i, i; 1; s), ^F,{l,hhl 



2-rii,g, 2) 3) J- — s) are solutions holomorphic in the neighborhoods of 
s = and s = 1 respectively, of the hypergeometric equation 



X 



3(s-r 



X 



1/12 
s(s - 1 



:X 



(33) 



with parameters a 



i b 

6' 



I and c 



1. The constants A,B are determined as before 



by considering the analytic continuation of the solutions near s = 0, where 



log(s)2Fi(|,i;l; 



n=0 



with hn '■= 2^/^(1 + n) — ^/'(| + ra) — '?/'(l + n). Then from fl32l) . one has 



= Ai(-log(s) + + 5, Ai = A 



r(i)r(i) 



where h{s) is a holomorphic function near s = 0. As indicated earlier, since z ^ z + 2 
stabilizes the vertex at z{0) = ioo, the pair {xi, X2} must form a basis of solutions near s = 
with monodromy (0 i)- Thus, the constant A is determined by taking Ai = i/ir above, 
so that z z + 2 onto the next branch as s makes a circuit around s = 0. The constant B 
is then obtained from the condition z{l) = e^'^*/'^, which yields B = — e'^*/^ after a similar 
calculation as in case 1. With these values of A, B in ( |32l) . it is possible to construct the 
power series for q := e'^*^'-^^ near s = 0, and its inverse s{z) = B2q{l + hiq + h2q^ + . . .) 
which gives the g-series for the triangle function S'(0, |, |; 2). The constant B2 = i'iS^/S, 
which follows from the leading coefficient = 2-?/'(l) — V'(l) ~ V'(l) above and evaluation 
of Digamma functions at rational arguments. 

From the Wronskian W{xi, X2) of the two solutions of 
C = Ai = i/iT in ([27]). With this value of C, and bi = 62 
gives the parametrization 



one computes the constant 
I from Table 1, ([27D then 



y{z{s)) = niP = Tci{l - s)3(xi + Gsxix'i) , Xi(s) 



2-^l(|' |i ^) 



One also deduces from arguments aligned with Proposition H] that. 



(34) 



Q{q) = {l-s)sxt, R{q) = (1 - ls)xl 



4.3. Cases 3, 4 & 6 



The automorphism groups in these cases correspond to the level- congruence subgroups 
ro(A^) of the modular group r(l), for = 2, 3, 4. The congruence subgroup ro(A^) is 
defined by 

To{N) := |7 = (^^ G SL2(Z) : c = (modN) 

The fundamental triangle T in each of these 3 cases has 2 parabolic vertices at z{0) = 
ioo, z{oo) = 0, while the vertex z{l) corresponds to an elliptic fixed point of order 2 at 
z(l) = i for ro(2), an elliptic fixed point of order 3 at z(l) = e^'^'^l'^ for ro(3), and a 
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parabolic point at z{l) = | for ro(4). Recall from Section 3 that in Case 6(b), the auto- 
morphism group is r(2) which is conjugate to ro(4). The principal congruence subgroup 
r(2) is defined as 

r(2) := |7 G PSL2(Z) : 7 = (J (mod2)| , 

which is a normal subgroup of r(l) of index 6. The fundamental triangle of r(2) is same 
as that of ro(4), but with respect to the variable z/2. That is, the parabolic vertices are 
located at z{0) = zoo, z{oo) = 0, and z{l) = 1 with respect to the local coordinate z eH. 

The triangle function S{0, P,0; z) in each of these cases will be constructed by de- 
veloping a g-expansion near the parabolic vertex z{0) = too, which is stabilized by the 
translation z ^ z + 1. From fl26al) . the hypergeometric equation corresponding to expo- 
nent differences 0, /3, is 

^ +U + 7^>+47(I^^ = °' ^^'^ 

with parameters a = b = (1 — /3)/2, and c=l. So the exponents are 0, at the regular 
singular point s = 0, and a, a at s = oo. A one-dimensional space of solutions for 
f l35|) near s = and s = oo respectively, are spanned by the functions 2-Fi(a,a; l;s) and 
(— s)~" 2Fi{a, a; 1; s"^). The second linearly independent solution near each singular point 
contains logarithms. The conformal mapping is defined as 

z{s) = A— , xi = 2Fi{a,a]l;s) , X2 = {-s) 2Fi{a,a]l; s ), (36) 
Xi 

where the constant A is determined by considering the analytic continuation of X2 near 
s = 0. This follows from the formula (see e.g. [35]). 



^ / In smlvra I , 

TT 



°° ( \2 

logi-shF,{a,a;l;s) + J2rhhns'' 



n=0 ^ ' 



— TT < arg(— s) < vr, and /i„ \=2'ip{l + n)—'tp{a + n)—%lj{l — a — n). Proceeding analogously 
as in the previous two cases, z{s) in (!32|) takes the form 

f \ w 1 / N 7 / A A sinCvra) 

z{s) = Ai{-\og{-s) + h{s)), Ai=A 



vr 



and h{s) is holomorphic near s = 0. The constant A (see Table 2) is determined by 
taking Ai = i/2'K above for cases 3, 4, and 6a, such that {xi, X2} has monodromy (01) 
at s = 0. Finally, inverting the power series for q := e^'^*^'^**-' leads to the expansion 
s{z) = Bq{l + hiq + h2q^ + . . .) for 5(0, /5, 0; z). The value of the constant B = —e^^° with 
ho = 2iIj{1) — il){a) — ■?/'(l — a), is listed in Table 2 below for each group ro(A^). For case 
6b corresponding to the group r(2), the parabolic vertex 2;(0) = ioo is stabilized by the 
translation z — > 2; + 2. Thus one needs to choose Ai =1/11 above, which leads io A = i 
and q = e^^^^^^ in this case. 

We remark that the triangle functions 5(0, /3, 0; z) associated with ro(A^), = 2, 3, 4 
as well as r(2) can be expressed compactly in terms of Dedekind eta functions. In fact. 
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N 


/3 


a=(l-/3)/2 


A = i csc(7ra)/2 




2 


1 

2 


1 

4 




-64 


3 


1 
3 


1 
3 


^/\/3 


-27 


4 





1 
2 


V2 


-16 



Table 2: parameters for constructing the triangle functions S{0, P,0; z). 

such expressions are available for the automorphic functions associated with all subgroups 
ro(A^) which are of genus zero, see Ref [26] for a complete list. 

For the cases 3, 4 and 6a, the Wronskian of the two solutions in (136|) is given by 
W{xi, X2) = Cs-\s - 1)-^", where C = (-l)2«/27ri. Then from ([27D one obtains the 
parametrization, 

y{z{s)) = -Qm{l - s)'"'-'[{{l + 2a - 61 - 62)5 + (62 - l)}xl + - , (37) 

where Xii^) = 2-^1 (a, a; 1; and (61,62) is given in Table 1. Note that in case 6b, the 
constant C = i/n in the Wronskian function. The parametrization of the Ramanujan 
functions P, Q, R for all three cases are listed in Table 3. These are obtained in the same 
manner as outlined in Proposition HI Note that there are two distinct parametrizations for 
each of the three cases corresponding to ro(A^), = 2,3,4. For ro(2) or ro(3), the pair 
of parametrizations are related via the involution s —>■ as noted below Table 1. Under 
this involution, it follows from fl5Bl) that z ^ ^ as the two linearly independent solutions 
Xi and X2 of fl5Sl) are switched in the quotient z. This is related (under appropriate 
normalization) to the Fricke involution: z —1/Nz for ro(A^) (see [26| for details). 
Thus in Table 3, the two parametrization in the case of ro(2) or ro(3) are transformed 
from one to the other by switching x = 2-^1(0-, a; 1; s) to x = 2-^1(0, a; 1; s) , s ^ 

as well as, taking into account the sign reversal in the Wronskian W appearing in (1271) 
through the chain rule formula d/dz = s'{z)d/ds = {x^ /W)d/ ds. 



r 


(61,62) 


P = yj-ni 


Q 


R 


ro(2) 


V2' 3/ 


2(\-sfl\x'\^sxx!) 


(4 - s)x' 


(1 -s)i/2(s + 8)x^ 


ro(2) 




(1 - sfl\x' + 125XX') 


(l-4s)x^ 


(l-s)i/2(8s + l)x' 


ro(3) 


(3' 2) 


3(s-l)^/^(x^ + 4sxx') 


(,_l)l/3(,_9)^4 


(27- 18s-s2)x^ 


ro(3) 


(14) 


{s - \fr\x' + I2sxx') 


(s - l)i/3(9s - l)x^ 


(1 + 18s -27s2)x^ 


ro(4) 


(14) 


(1 - 2s)x^ - 12s(s - \)xx! 


16(s^ - s + ^)x^ 


32(2s-l)(s^-s-4)x^ 


r(2) 


(14) 


(1 - 2s)x^ - 6s(s - \)xx! 


{s^-s + l))x' 


i(2s-l)(s^-s-2)x^ 



Table 3: Parametrization of P, Q, R. Here x = 2-^1(0, cl', 1; s) and s{z) = S{0, (3, 0; z) with 
a and /? values given in Table 2. 

The relationship between the parametrization of the Ramanujan functions (Eisenstein 
series) P, Q, R listed in Table 3 and Ramanujan's theories of elliptic integrals will now be 
established. As mentioned in Section 1, Ramanujan originally gave the parametrization 
d?]) in terms of complete elliptic integrals of the first kind. These correspond to the 
parametrization associated to the groups ro(4) and r(2) in Table 3. The parametrization 
corresponding to ro(2) and ro(3) are related to the parametrization ([8]) in Ramanujan's 
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alternative theories. The hypergeometric functions appearing in ([7]) and ([8]) are related 
to those in Table 3 via the well-known Pfaff transformation 



2Fi{a,b;c; s) = {1 - s) " 2Fi{a,c - b;c;x{s)) , x{s) 



s - 1 

so that in the present case with a = b, c = 1, (136|) takes the form 

2Fi(a,a;l;s) 2sm(7ra) 2-ri(a, 1 — a; 1; x) 

Then by setting a = 1/r and 2TTiz{s) = —Ur, ( !38l) coincides with ([8]), and thus one 
recovers from Table 3 the Eisenstein series parametrizations in Ramanujan's theories for 
signatures r = 4,3,2 corresponding to the groups ro(2),ro(3) and ro(4), respectively. 
For completeness, these are presented in Table 4 below, in terms of the hypergeometric 
function Xr{x) := 2-^1(^5 which appears in Ramanujan's theories. 

4.4- Case 5 

The triangle group corresponding toa = 0,/3 = |,7 = 0is not a Fuchsian group of first 
kind is not a positive integer) although the triangle functions s{z) = 5(0, |, |; z) and 

's{z) = S{0,^,0;z) are related via s{ez) = —A's{z)/{1 — 5(2;))^, e = —\J\^ as mentioned 
in Section 3. The main significance of this case lies in the fact that it corresponds to 
Ramanujan's alternative theory of signature r = 6 as outlined below. 

The triangle function S'(0, |,0;z) can be constructed employing the hypergeometric 
theory in exactly the same way as for the cases associated with the subgroups ro(A^), = 
2,3,4 of the modular group, even though this case is not known to have any "modular 
interpretation" [26]. Specifically, one proceeds from the hypergeometric equation (!35l) 
with /5 = |, and defines z{s) as in (l36l) with a = i. Then one finds that the constants 
A = i and B = —e^° = —432. The Wronskian in this case turns out to be W{xi,X2) = 
Cs-\s - l)-i/3 ^i^g^g Q ^ -l/2TTi. Finally, substituting (61,62) = (|, |) ^^"^ ^^^^^ ^ 
into (!37|l and using arguments similar to that in Proposition HI one obtains the following 
parametrization for the Eisenstein series in Ramanujan's theory of signature r = = 6, 



y{z{s)) 1 . 2 

ni 

1 2Fi(l,i;l;.-i 



P{q) = = (1 - s)Hxl + I2sxx') , Qiq) = (1 - s)sx' , Riq) = (1 + s)x' 



Xis) = 2F,il i; 1; .) q = exp - 27r(-.)-6 / J . (39) 

The map in ( l38l) induced by the Pfaff transformation also applies to this case, and leads 
to the alternative parametrization for r = 6 listed in Table 4. 

The parametrizations listed in Table 4 can be found elsewhere, e.g., in Refs. [2^ [25] . 
The first entry in Table 4 for r = 2 appears in ([7j) as Ramanujan's original parametrization. 
The second entry (case 2*) corresponds to the parametrization of the Eisenstein series in 
terms of the elliptic modular function x = X{z). Recall from Table 1 that the triangle 
function for r(2) is S'(0, 0, 0; 2;/2) where S'(0,0,0;2;) is the triangle function for ro(4). 
Then from the formula given at the end of Section 3, it follows that S'(0, 0, 0; 2;/2) = 
X[z)[X{z) — The expression of x{s) in the Pfaff transformation formula implies that 
the involution s —>■ - corresponds to the involution x — * 1 — x, as well as the switching of 
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r 


(&1,&2) 


P = yj-Ki 


Q 


R 


2 


f- -) 


(1 - 5X)X2 + 12x(l - X)X2X2 


(1 + 14x + 


(1 + a;)(l - Ux + x'^)xl 


2* 


M- 


(1 - 2x)xl + 6a:(l - x)x2X2 


(1 + 


{l + x){l-^x + x')xl 


3 




(3 - Ax)xl + 12x(l - x)x3X3 


(9 - 8x)xt 


(8x2 _ ^ 27)xl 


3 




(1 - Ax)xl + 12x(l - x)x3X3 


(1 + 8x)x^ 


(1 - 20x - 8x'^)xl 


4 




(2 - 3x)xl + I2x{l - x)xix'A 


(4 - ?,x)x\ 


(8 - 9x)x^ 


4 




(1 - 3x)x4 + 12x(l - x)x4Xl 


(1 + ?>x)xl 


(l-9x)x^ 


6 




(1 - 2x)xl + 12x(l - x)x6X6 


xt 


(1 - 2x)x^ 



Table 4: Parametrization of P,Q,R via Xr = 2-Fi(-, 1; a;), = s/{s — 1). The 
corresponding triangle functions are s{z) = S{0,P,0; z) with /3 = 1 — ^. 



the two 2F1 functions in the quotient z{s) in (138|) . This leads to the similar transformation 
as in Table 3, between the two cases for r = 3 or 4 listed in Table 4. It is interesting 
to note that the parametrization in the sextic case remain invariant under the involution 
X ^ 1 — X. 

5. Concluding remarks 

In this note, we have reviewed the relationship between the Chazy equation and the con- 
formal maps defined by ratios of solutions of Fuchsian equations with 3 regular singular 
points. The Chazy equation turns out to be a particular case of more general third or- 
der nonlinear differential equations which can be derived systematically by developing 
the transformation properties of functions under the projective monodromy group asso- 
ciated with the Fuchsian equation. Much of this note has been devoted to exemplify 
the important connection between the Chazy equation and Ramanujan's study of elliptic 
integrals and theta functions which play a defining role in the contemporary theory of 
modular forms and elliptic surfaces. By the way of elucidating this beautiful relationship, 
we have derived all possible parametrizations of the Chazy solution via the Schwarz 
triangle functions. We show that these parametrizations are also related to Ramanujan's 
parametrization for the Eisenstein series P,Q,R in terms of hypergeometric functions, 
arising in his theories of modular equations of signatures 2, 3, 4 and 6. We also give two 
additional hypergeometric parametrizations stemming from the triangle functions associ- 
ated with the modular group r(l) and its index 2 subgroup F^. The case corresponding 
to the full modular group F(l) was considered by Chazy himself Furthermore, we 
note that it is possible to systematically construct a number of third order nonlinear 
equations associated with modular as well as other automorphic groups. Some of these 
equations corresponding to the subgroups groups Fo(2), Fo(3) and Fo(4) of F(l) have been 
studied recently in Refs. El] and in Ref. which analyzes the Ramanujan's differ- 
ential system (jlj) as mentioned in Section 1. We end this note with an application of the 
parametrization of Eisenstein series in Table 4 to solutions of differential equations. 

Consider the DH system of differential equations which was introduced in Section 
2, and which corresponds to the system (IT^ with r = 0. As mentioned there, if the 
variables Wj, j = 1,2,3 satisfy the DH system then y{z) := 2{wi + W2 + W3) solves the 
Chazy equation ([2]). In the process of verifying this assertion, one computes y'iz) = 
2{wiW2 + W2WS + WsWi) and y"{z) = 12wiW2Ws, which implies that the DH variables Wj 
are the distinct roots of the cubic — {y/2)w'^+ {y'/2)w — y"/l2 = 0. Expressing y, y', y" 
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in terms of the functions P,Q,R (see e.g. Proposition Hj) and introducing t = {Q/m)w, 
one obtains the cubic equation 

- SPt^ + 3(P - Q2)t - (P^ - 3PQ + 2P) = , 

whose coefficients are polynomials in P, Q, R- Now note from Table 1 that y{z) := 2(toi + 
W2 + wz) corresponds to the case 6b associated with the group r(2). The hypergeometric 
parametrization for this is the case r=2* in Table 4. Using those parametrizations for 
P, Q, R in the above cubic, one can verify that the cubic factorizes as {t — ti){t — t2){t — t^) 
with 

ti = 3(1 - a;)x2 + ^3 , t2 = -?>xxl + h, h = Qx{l - x)x2X2 ■ 

Thus one immediately recovers the special solutions Wj = {ni/6)tj, j = 1, 2, 3 for the DH 
system in terms of elliptic modular form x = X{z) and X2 = 2Fi{^, |; 1; x) = {2/7t)K{x), 
where K{x) is the complete integral of the first kind defined in Section 1. Such lin- 
earization procedures to derive exact solutions can be also employed to a large class of 
nonlinear differential equations with automorphic properties as discussed in this article. 
These problems will be addressed elsewhere. 
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